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The main goal of the present paper is to study how polarization of photons affects their motion 
in a gravitational field created by a rotating massive compact object. We study propagation of 
the circularly polarized beams of light in a stationary gravitational field. We use (3+l)-form of 
the Maxwell equations to derive a master equation for the propagation of monochromatic electro- 
magnetic waves of the frequency uj with a given helicity. We first analize its solutions in the high 
frequency approximation using the 'standard' geometrical optics approach. After that we demon- 
strate how this 'standard' approach can be modified in order to include the effect of the helicity 
of photons on their motion. Such an improved method reproduces the standard results of the ge- 
ometrical optics at short distances. However, it modifies the asymptotic behavior of the circularly 
polarized beams in the late-time regime. We demonstrate that the corresponding equations for the 
circularly polarized beam can be effectively obtained by modification of the background geometry by 
including a small factor proportional to ui^ 1 whose sign corresponds to photon helicity. We obtain 
the modified equations for circularly polarized rays by using such a 'renormalization' procedure, and 
calculate the corresponding renormalization term for the Kerr geometry. 
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I. INTRODUCTION 

The main goal of the present paper is to study how 
polarization of photons affects their motion in a gravita- 
tional field created by a rotating massive compact object. 

It is well known (see e.g. [l|-y]) that the electromag- 
netic field equations in an external gravitational field are 
formally equivalent to the Maxwell equations in a flat 
spacetime written in Cartesian coordinates in the pres- 
ence of a "medium" with the dielectric permittivity and 
magnetic permeability tensors related to the spacetime 
metricQ. This approach was used to study different as- 
pects of propagation of electromagnetic waves in a grav- 
itational background d 04^1- I n particular, using this 
approach it was shown that for the motion of a photon 
in the external stationary gravitational field its helicity is 
conserved dUl. This implies that in a stationary exter- 
nal gravitational field, in the absence of photon creation, 
the state of the circular polarization remains the same. 
In particular, scattering of the incoming right (left) cir- 
cularly polarized electromagnetic radiation in an asymp- 
totically flat spacetime results in the outgoing right (left) 
circularly polarized radiation. For such problems one 
can decompose an electromagnetic wave into two inde- 
pendent non-interacting components of the left and right 
circular polarization^. 



* vfrolov@ualberta.ca 
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1 It is interesting to notice that the Maxwell equations in a curved 
spacetime filled with a moving fluid with non-trivial dielectric 
permittivity and magnetic permeability can be identically rewrit- 
ten as the vacuum Maxwell equations in a specially modified 
so-called Gordon's optical metric. This result was obtained by 
Gordon 4]. For recent discussion and applications see, e.g., Q 
and [fj. 

2 If a gravitational field creates photons, an initially circularly po- 



If the wavelength of electromagnetic waves is much 
smaller than the characteristic scale of the gravitational 
field, one can use the geometrical optics approximation 
for their description. In this approximation the wave 
equations are reduced to the Hamilton- Jacobi equation 
for the eikonal and transport equations for the amplitude 
of the wave. The former is equivalent to geodesic equa- 
tions for null rays. If an electromagnetic wave is linearly 
polarized, then in the geometrical optics approximation 
the vector of polarization isparallelly propagated along 
the ray (see, e.g., and [ll|, p. 570). First order cor- 
rections to the geometrical optics approximation where 
discussed by Ehlers [l(| and polarization effects in the 
Schwarzschild metric were studied in [l2| and [l3j . 

In a general result of the bending of a light 

ray, the polarization vector changes its direction. Be- 
sides this, if the source of the gravitational field is rotat- 
ing, the polarization vector can rotate around the prop- 
agation vector. This effect is a gravitational analogue 
of the well-known electromagnetic Faraday effect, which 
appears in magneto-active media [3 [l5[ . This effect was 
studied in the Kerr spacetime, where the equations of the 
parallel transport along null geodesies can be solved ex- 
actly [l6|-Qm . The relativistic effect of Faraday rotation 
may be important for the explanation of the properties 
of radiation from accreation disk of a rotating black hole. 
In particular, Laor, Netzer, and Piran demonstrated that 
this effect can suppress polarization of the optical radia- 
tion from the active galactic nuclei [201 ] . 



larizcd beam of photons may have in the final state created pho- 
tons of different polarization. This can happen if the gravita- 
tional field is time dependent, or stationary, for example in the 
case when the beam is propagating through the ergosphcrc of a 
Kerr black hole, where the superradiant paticle creation is pos- 
sible. 
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The effect of the gravitational Faraday rotation of the 
spin of a photon is the result of the action of the pho- 
ton's trajectory in a given gravitomagnetic field on the 
spin orientation. One can expect a dual effect, namely an 
action of the helicity of the photon on its trajectory. Such 
an effect is analogous to the optical Magnus effect which 
was studied for the propagation of light in an optically 
inhomogeneous medium [21] . In other words, when a po- 
larized radiation propagates near a rotating gravitating 
body photons of the opposite initial helicities are scat- 
tered differently [1, [22[ . It was also shown that the scat- 
tering amplitude of the left and right circularly polarized 
waves are different 0, Because of these effects, one 
can expect a partial polarization of the scattered by the 
Kerr black hole initially unpolarized radiation Q. The 
helicity asymmetry due to the gravitational deflection of 
unpolarized light by a rotating body was calculated in 

dim. 

The dependence of photons scattering by a rotating 
black hole on their helicity is an effect beyond the lowest 
order of the geometrical optics approximation. Suppose 
that two beams of light of the same wavelength A but 
opposite helicity pass near a rotating object of the mass 
M and the angular momentum J at the distance L 3> 
GM/c 2 . Then the expected angular difference of their 
asymptotic trajectories is of the order of GJX/(c 3 L 3 ). 
This estimation was obtained by Mashhoon [t| [22| in the 
weak field approximation for the gravitational field of a 
rotating body. 

The aim of this paper is to study the propagation of 
circularly polarized beams of light in a stationary space- 
time assuming that the gravitational field is not weak. 
This paper is organized as follows. In Sec. II we dis- 
cuss general aspects of the geometrical optics approxi- 
mation for the polarized light propagation in an external 
stationary gravitational field. This section contains ref- 
erences to the numerious results related to our problem 
and its analogues for the spin particle motion in a mag- 
netic field. In Sec. Ill we use (3+l)-decomposition of 
a stationary gravitational field and derive 3-dimensional 
master equation for the propagation of monochromatic 
elecromagnetic waves with a given state of helicity. The 
'standard' scheme of the geometrical optics approxima- 
tion applied to the master equation is presented in Sec. 
IV. Section V describes a modified approximation which 
takes into account the back reaction of the photon polar- 
ization state on its trajectory. We also demonstrate that 
the corresponding modification of the photon equation of 
motion can be effectively obtained by simple 'renormal- 
ization' of the background gravitational field, depending 
on the polarization of the photon and its frequency. In 
order to illustrate this procedure, we calculate the corre- 
sponding 'renormalization' term for the Kerr geometry. 
General remarks and possible applications of the devel- 
oped formalism are discussed in Sec. VI. The Appendix 
contains additional technical results. 

In this paper we use the system of units where G = 
c = 1 and the sign conventions for the metric and other 



geometrical quantities adopted in the book by Misner, 
Thorne, and Wheeler [ll|. 



II. APPLICATION OF THE GEOMETRICAL 
OPTICS TO THE POLARIZED LIGHT BEAM 
PROPAGATION 

Description of the light propagation by means of the 
geometrical optics is an approximation. This approxi- 
mation is valid when the wavelength of light A is much 
smaller that other characteristic scales of the problem, 
such as (i) the length l m at which the parameters of the 
medium change significantly; (ii) the radius Ir of the 
spacetime curvature, and (iii) the radius l w of the curva- 
ture of the wave front surface. Let I be the smallest of 
these scales, then the condition of validity of the geomet- 
rical optics approximations is e = X/(2nl) <C 1. Under 
this condition one expects that the field is fast oscillating 
while its amplitude changes slowly. The mathematical 
version of these physical ideas are asymptotic formulas 
for solutions of the corresponding field equations. The 
method of obtaining these asymptotics is well known (un- 
der different names) and widely used both in physics and 
mathematics. For example, in quantum mechanics it is 
called the WKB method, and the smallness parameter e 
is proportional to the Planck constant h. 

The general scheme of the construction of the asymp- 
totic stationary solutions in the geometrical optics ap- 
proximations is the following (see, e.g., [HI): 

• Let F(x) be a field defined on an n-dimensional 
manifold with coordinates x. Write it in the 
form F{x) = f(x) exp(iS (x) / e) , and substitute this 
ansatz into the field equation(s); 

• In the obtained relation select the leading in e 
terms. A condition that the corresponding expres- 
sion vanishes gives the eikonal equation, which is 
usually of the form if(V«S, a;) = 0; 

• In order to solve this first order partial differen- 
tial equation, one considers it as the Hamilton- 
Jacobi equation for a particle with the Hamiltonian 
H(p,x), where p = V5 is the particle momentum, 
and the initial conditions S = Sq{x) — const., such 
that po(x) = VSq(x). 

• Let £ be a parameter along the particle trajectory, 
and let x(£)\e=o = xq, p(£)\e=o = po be the initial 
conditions. 

• Solve the Hamilton equations of motion for the ini- 
tial conditions. The set of the solutions {x(£),p(£)} 
determines an n-dimensional submanifold of the 
phase space, which is known as a Lagrange man- 
ifold. Its characteristic property is that the sym- 
plectic form vanishes on it. The function S(x) is 
the action defined on the particle trajectories. 
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• Next to the leading order in e equations give trans- 
port equations which allow one to find the ampli- 
tude F along the particle trajectories, and hence, 
determine it on the Lagrange manifold. At the 
points where the projection of the latter on the 
configuration space is regular, this projection de- 
fines the required asymptotic expansion of the am- 
plitude f(x) ~ / (a;) +efi(x) + .... 

• Points of the Lagrange manifold where its projec- 
tion on the configuration space is singular form 
caustics, which require special consideration (for 
details see, e.g., [26|,[27} an d the references therein). 

In the application of the geometrical optics to the prob- 
lem of the polarized light beam propagation several new 
features become important. First of all, the light ampli- 
tude has not one but several components. As a result, the 
transport equation for the amplitude propagation along 
the particle trajectories has a matrix form. One of the 
methods used in such case is based on the diagonaliza- 
tion of the corresponding equations (see, e.g., |28l. |29|). 
A similar approach is used for the construction of the 
geometrical optics asymptotic solutions of the Maxwell 
equations in magneto-active media. In this case, after the 
diagonalization one has independent modes. Moreover, 
the eikonal equations for different modes are generally 
different (see, e.g., [3(|[3l[ and the references therein). A 
remarkable desription of the general eikonal method for 
systems with several components is given in [32J . 

Our main subject, spinoptics in a stationary gravita- 
tional field, in many aspects is similar to more analyzed 
and better understood problem of a charged spin parti- 
cle motion in an external electromagnetic field. A well- 
known example of such a problem is the Stern- Gerlach 
experiment. Observations show that when a quantum 
charged particle with spin (electron) propagates through 
a region of inhomogeneous magnetic field it gets deflected 
differently for different spin orientations. In order to de- 
scribe this effect, one can use the WKB approximation 
to derive quasi-clasical equations of the electron motion. 
The problem is that the electron magnetic moment inter- 
acting with the gradient of the magnetic field is propor- 
tional to H, and hence, it vanishes in the limit h — > 0. As 
a result, a simple application of the WKB does not allow 
to obtain the required result for the spin-dependent de- 
flection of the electron in the Stern-Gerlach experiment. 
This conclusion made by Paul i f33j and supported by 
more accurate calculations (34l - l36j ] stimulated many in- 
teresting discussions (see, e.g., the review [3(| and the 
references therein) . 

A general resolution of this puzzle is the following. In 
the Stern-Gerlach experiment one deals with a scattering 
problem. That is why in the quasi-classical description 
one needs to know the WKB asymptotic solution not 
only locally, but globally as well. The spin-dependent 
terms, which are small (of the order of K), locally give a 
small contribution to the electron trajectory. But as a 
result of the accumulation of these spin-dependent cor- 



rections during the scattering, their effect on the electron 
trajectory at the later time is not small at all. In order 
to be able to describe this effect and to construct the 
corresponding global WKB solution one needs to include 
the lowest order spin-dependent corrections in the eikonal 
equation. This results in the global modification of the 
corresponding Lagrange manifold. One can expect that 
this modification is locally small, but becomes important 
in the late-time regime. This approach was originally ad- 
vocated by de Broglie [37} ■ It was made more concrete 
ten years later by Schiller [H, . 

Several different approaches to the problem of motion 
of a spin 1/2 particle in a magnetic field in the WKB 
approximation have been developed recently. The most 
consistent approach is to include spin degrees of freedom 
into the quasi-classical equations (see, e.g., (40l - |42l |). As 
a result, the phase space of the particle becomes a di- 
rect product the standard 2n-dimensional phase space 
and a two-dimensional sphere. The Hamiltonian defined 
on this enlarged phase space describes both the spin ro- 
tation induced by the electron motion in the magnetic 
field and the back reaction of the spin rotation on the 
particle motion 42] . In this scheme, the symplectic form 
is also modified and is determined on the enlarged phase 
space 43]. A similar geometrical approach to a theory 
of spinoptics in refractive media can be found in [44| and 

53- 

Another way to include spin-field interaction is to mod- 
ify the eikonal by including into it a term dependent on 
the spin orientation, which is proportional to h [46] . This 
modification does not change the local WKB expansion of 
the field and can be considered as its trivial redefinition. 
However, these locally small corrections of the eikonal 
equation result in the non-trivial change of the late-time 
asymptotics of the particle trajectories. A similar result 
was obtained in |47| where there was formulated a general 
principle which can be used for identification of the terms 
which must be included into the action S. The authors 
called this method the modified geometrical optics. 

An analogical effect of the photon spin on its motion is 
known as the spin- Hall effect [48M52I ] . This intrinsic spin- 
Hall effect, associated with the Berry phase, was observed 
experimentally for the propagation of helical light beam 
at a grazing angle inside a glass cylinder [48| . The Hamil- 
tonian and Lagrangian approaches to the problem of evo- 
lution of spinning light in an inhomogeneous medium are 
summarized, respectively, in the papers [53| and [sij . 

In our approach to the spinoptics in a stationary grav- 
itational field, we demonstrate how the next to the lead- 
ing order helicity-dependent corrections naturally arise 
in the eikonal equations and develop the corresponding 
improved scheme, which allows one to describe the ef- 
fect of the state of circular polarization of the light beam 
on its motion. We use the similarity of this problem to 
the better-understood case of the spin particle motion 
in a magnetic field. In spite of this similarity, the case 
of light motion has one essential difference: photons are 
massless while electrons have mass. In the quantum de- 
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scription of a photon, its state can be characterized by 
the photon's helicity. The spin vector of massive par- 
ticles can have an arbitrary orientation with respect to 
their velocity vector. In the classical description of light, 
a state with a given helicity corresponds to either right 
of left circular polarization of the electromagnetic wave. 
The photon's helicity is a conserved quantity 0,0. If 
an external gravitational field does not create photons, 
the number of photons in a beam and their helicity re- 
main the same in the scattering process. In other words, 
solutions of the source-free Maxwell equations in an ex- 
ternal gravitational field can be separated into the sum 
of two solutions, one with right and another with left cir- 
cular polarization. In a general case, this decomposition 
involves self-dual and anti-self-dual complex fields (55[. 
Thus, to study the spinoptics effects in a gravitational 
field, it is natural to consider electromagnetic waves with 
a given circular polarization. In our approach, we write 
such decoupled wave equations and apply the geometri- 
cal optics approximation to each of the two sectors corre- 
sponding to the right and left circularly polarized waves. 

However, there is some problem remains. Namely, the 
Maxwell equations form a set of eight first order par- 
tial differential equations for six field variables, which are 
components of the electromagnetic field tensor . This 
redundancy is the consequence of the gauge invariance 
of electromagnetic field. Using (3+l)-decomposition we 
write the complete set of the Maxwell equations in an 
equivalent form of the master equation. The master 
equation follows from the set of three first order partial 
differential equations for 3-dimensional complex vectors 
of the field components with a given state of polariza- 
tion. The form of the master equation is used to develop 
the scheme of the geometrical optics approximation. We 
shall also demonstrate that the scheme is convenient and 
provides a natural starting point for the new improved 
scheme, which is the purpose of our work. 



III. PROPAGATION OF CIRCULARLY 
POLARIZED LIGHT IN A STATIONARY 
SPACETIME 

In this section we consider the propagation of a weak 
electromagnetic field in a stationary gravitational field 
background. In order to fix our notations we collect the 
required equations. 



A. (3+l)-decomposition of a stationary metric 



lines is arbitrary. Thus, there exists an arbitrariness in 
the choice of t: 



t-tt = t+ q{x { ) , 



(1) 



where q{x l ) is some differentiable function of the coor- 
dinates x 1 . One can choose the Killing time t and the 
coordinates x l labeling the Killing trajectories as space- 
time coordinates (t = x° ,x l ). 

In these coordinates the spacetime metric can be pre- 
sented in the following form: 



ds 2 



h(dt - gi dx 1 ) 2 + hjij dx % dx 3 , (2) 



where h = — £ 2 > 0. In what follows, all the opera- 
tions involving the Latin indices are performed 
by using the metric 7y and its inverse j lJ . In particu- 
lar, the covariant derivative in this metric is denoted as 
Vj(...) = We preserve the notation (. . .) ;jU for the 
4-dimensional covariant derivative defined with respect 
to the metric g^,. 

The Killing equation + £„ ;ft = implies that the 
metric functions in are time- independent. Under the 
coordinate transformation (fT]) h and Ty remain invariant, 
while gi's transform as 



9i 9i 



9i 



(3) 



It is possible to use the gauge transformation (fTJ to put 
g t . i = 0. To find this transformation it is sufficient to use 
q(x l ), which is a solution of the equation q'- l A = —g 1 . i - 

The metric tensor 7 M „, which in a choosen coordi- 
nate system has the components 7^,, = S^S^jij, is pro- 
portional to the projector onto a surface orthogonal to 



0- 



e 



(4) 



In a static spacetime, i.e., when gi — 0, the metric Tj, 
reduces to the so-called optical metric (see, e.g., [56H58| 
and the references therein) . Spatial trajectories of light 
rays are geodesies of the optical metric. 

The metric functions in the (3+l)-decomposition of 
the metric ([2]) are related to the 4-dimensional metric 
functions g^ v as follows: 

h = -goo , 9i = -j^ , -9 = ^ 4 7 j (5) 



A stationary spacetime possesses timelike Killing vec- 
tor generating the time-symmetry transformation. 
The integral lines of this vector field obey the equation 
dx^/dt = We use the Killing time t as a parameter 
along the integral lines of £ M and introduce three other 
coordinates x l , i — 1,2,3, which are constant along the 
lines. The origin t = of the Killing time on each of the 



9ij 9oi9oj _ 9ij 



9i9j ■ 



(6) 



5oo #00 ^ 

Here g = det(g M „) and 7 = det(7y). Let us denote g % = 
^9j- Then using the relation g^g = 6* one obtains 

7«=V 3 ", 9 l = hg°\ g 00 = -h- 1 (l-9i9 i ). (7) 
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B. Maxwell equations in (3+l)-form in a 
stationary spacetime 

Consider a test electromagnetic field in a stationary 
spacetime. In accordance to the (3+l)-decomposition of 
the metric ([2]) one can also make a similar decomposition 
for the components of the electromagnetic field. Follow- 
ing [H| we denote 



Ei 



1 i0 : 



R = F 



(8) 



D i = h 2p0i ^ jjij = tfpij , (9) 
These objects are not independent. Using the relations 
Fio = g^gouF^ , F* = g^g^F^ (10) 

we derive 

Di = Ei- Hijg j , I!" //'•' f.y • E j g\ (11) 

Let Cijk and e 1 ^ be 3-dimensional completely antisym- 
metric Levi-Civita symbols normalized by the conditions 
£123 = 1 and e 123 = 1, and 



e ijk = y/ie ijk , e l]k = -^-e ijk 
V7 



(12) 



be the corresponding completely antisymmetric tensors. 
These tensors obey the relations 



eijk-.l = , e jk :l = . 
We introduce the following vectors: 

B i = i, «* B jk , B l = X -e ldk W k 

Using the relation 



Imk zlxm X?ncl 



we obtain 



B l} = e ijk B k , //'•' = <"'• //,. . 



(13) 



(14) 



(15) 



(16) 



To define a vector product of vectors A and B in a 3- 
dimensional space we shall use the following notation: 



where 



For example, 



C = [A x B] , 

C l = e iih AjB k . 
1 



(17) 
(18) 



[H x g y = -e^ k e jlm B lm g k = B lk g k . (19) 

Thus the relations (fTTj) take the form 

D = E — [g x H], B = H + [g x E] . (20) 



We shall introduce the following operations defined in a 
3-dimensional space: 

1 



divA=(V,A) = A\ A = — 

V7 



(21) 

(curL4) 1 = [V x A} 1 = e ijk A k:j . (22) 
It is easy to check that the following relations are valid: 



div(M) = (V/,A) + /divA, (23) 

curl(/A) = [V/ x A] + /curLA , (24) 
curl [A x B] = AdivB - BdivA 

+ (B,V)A-(A,V)B, (25) 
V(A, B) = [A x curLB] + [B x curLA] 

+ (A,V)B + (B,V)A, (26) 

div[A x B] = (B, curLA) - (A, curlB) , (27) 

div curlA = . (28) 

To prove the last relation we notice that its left hand side 
can be written as 



e ijk WiWjA k = e ijk Ri jkl A l . 



(29) 



Since the Riemann curvature tensor obeys the identity 
R\ijk]i = 0, this expression vanishes. The relations (|23[1 - 
(|28p in a curved 3-dimensional space with the metric 7^ 
have the same form as the corresponding relations in a 
flat 3-dimensional space. 

Let us now present the source-free Maxwell equations 
in a (3 + l)-form. The first set of the Maxwell equations 



F[/j,u,X] — 
in a stationary spacetime gives 

divB = , curLE 



B . 



Here we denote A = Aj. 

The second set of the Maxwell equations 



1 



0, 



in a stationany spacetime gives 

divD = , curLff = D . 



(30) 
(31) 

(32) 
(33) 



The energy-momentum tensor of an electromagnetic field 



T, 



F F ' 



1 



-.g^v^x 



Act 



satisfies the relation 



T» v . v = . 



Thus, for the Killing vector field £^ = 5^ Q we have 
1 



--(V=g^Tn^ = 0- 



(34) 



(35) 



(36) 
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This relation implies 

£ + divV = , (37) 

where 

£ = ^[(E,D) + (B,H)], (38) 
is analogical to the electromagnetic energy density, and 
V=^[ExH], (39) 
is analogical to the Poynting vector. 

C. Master equations for circularly polirized light 

The system of equations (|3Tj) and (f3"3"|) can be rewritten 
in a complex form. For this purpose we introduce the 
following complex quantities: 

F ± = E±iH , G ± = D±iB. (40) 

In these notations the Maxwell equations (|3"Tj) and (f3"3"|). 
and the relations ([20")) take the form 

divG* = , curlF* = ±iG ± , (41) 

G ± = F ± ± i[g x F*] . (42) 

In a stationary space a monochromatic wave can be 
written as 

E = e~ luit e + e lut 6 , H = e~ iut U + e iut U , (43) 

D = e-^T) + e iclJt X> , B = e^B + e iu)t 3 , (44) 

where the bars over the vectors denote complex conju- 
gate. Substituting these relations into Eq. (|40p we derive 

F ± = e -<utjr± + e ^tj^ t (45) 

G ± = e -iutg± + jutgT _ ^ 

where 

3* =S±iH, g ± = T> ± iJ5 . (47) 

The complex positive frequency amplitudes and G^ 1 
describe left (for — ) and right (for +) circularly polarized 
waves. In other words, for the right (left) circularly po- 
larized monochromatic radiation the quantities with the 
label — (+) vanish. 

For the monochromatic waves the system of Eqs. I|4ip- 
reduces to 



curLF* = ±u J g ± , (49) 

g ± =J r± ±i[gxF ± }. (50) 

Equation (|4"8f is not independent. It follows from Eq. 
P9l> and the identity pg|). 

Equations (|49)l and (|50|) imply the "master equation" 

curLT^ = + ioj[g x J^] . (51) 

Solving this master equation for !F ± , we can find Q ± 
from Eq. (|49|) . and Eq. (l48l) is satisfied automatically. 
The form of the master equation implies that the left and 
right circularly polarized waves are decoupled. 

IV. 'STANDARD' GEOMETRICAL OPTICS 
APPROXIMATION 

A. Eikonal ansatz 

Let us briefly describe now the 'standard'geometrical 
optics in its 3-dimensional form is adopted to circularly 
polarized light propagation. Our main assumption is that 
the characteristic wavelength A = 2ttlo~ 1 is much smaller 
than any other characteristic length scale I of the prob- 
lem. We shall use cj _1 as a small parameter. In fact, it 
always enters the expressions in the dimensionless form 
e = (ojI)^ 1 . In what follows, we shall use w _1 as a small 
parameter in our geometrical optics expansion. Follow- 
ing to the geometrical optics construction (see, e.g., [Tpj 
and [U, p. 570), we split a monochromatic wave into a 
rapidly changing phase and slowly changing amplitude. 
Such a split corresponds to the following eikonal ansatz: 

T" = f a e lu)S . (52) 

Here a = ±1 defines polarization of the wave, as defined 
above, / CT is the complex valued amplitude, and S is the 
real valued phase. We shall always consider a monochro- 
matic wave of a fixed polarization, either a = +1 or 
(7 = — 1 , and skip for brevity the superscript cr, which we 
shall restore in the final expressions. 

Substituting the eikonal ansatz ([52")) into Eq. (T5"Tj) and 
using Eq. (|24[) we derive 

Lf = au}- 1 curl f . (53) 

Here L is a linear operator acting in a 3-dimensional com- 
plex linear vector space, 

Lf = f ia[n x /] , (54) 

where n = p — g and p = V5. The vector up has the 
meaning of the usual wave vector. 

According to the standard prescription, the slowly 
changing complex amplitude / can be expanded in in- 
verse powers of U) as follows: 



div£P = , 



(48) 



/ = /o+^ 1 /i+^ 2 / 2 + 



(55) 
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Here / is the leading order term in the geometrical op- 
tics approximation, and the higher order expansion coef- 
ficients /j's, i > 1 define post-geometrical optics correc- 
tions. Substituting the expansion (|55[) into Eq. (|53p we 
obtain the following equation: 

Lf + u~ 1 [Lf 1 - ercurl/ ] 

+ cj- 2 [Lf 2 - crcurl/i] + . . . = . (56) 

B. Properties of the operator L 

The leading order equation in the high frequency limit 

is 

Lf =0. (57) 

In an explicit form this equation reads 

L i j fi=[S l j -iae i kj n k ]f =Q. (58) 

In order to have a non-trivial solution / , determinant 
of L must vanish. The corresponding relation gives 

{n,n) = l. (59) 

This is the eikonal equation. It is easy to check that 
Ln = n. That is, n is an eigenvector of L with the 
eigenvalue 1. 

To determine other eigenvalues of L it is sufficient to 
find solutions of the characteristic equation 

det[(l - X)S i j - iae l kj n k ] = . (60) 

Calculating the determinant (see Appendix A) one ob- 
tains the following equation: 

(1 -A)[(l- A) 2 - (n,n)] =0. (61) 

Using Eq. ([59]) one finds that the other two eigenvalues 
are equal to and 2. 

To find the eigenvectors for A = 0, 2 let us consider two 
unit real vectors e.\ and e 2 such that 

(e 4 ,e i ) = l, (ei,e 2 ) = (e. n n) = 0, i = 1, 2 , (62) 
[ei x e 2 ] = n , [n x e{\ = e 2 , [n x e 2 ] = — e x . (63) 

The set of vectors {n, ei,e 2 } forms a right-oriented or- 
thonormal basis in a 3-dimensional linear vector space. 
Let us now define the following vectors: 

m= -^=(e 1 +i<je 2 ) , fh = -^=(e x - iae 2 ) . (64) 

The vector m is the complex conjugate of m. According 
to Eqs. (j6"2"j) and these vectors satisfy the following 
relations: 

(m, m) = (m, n) = , (m, m) = 1 , (65) 
[m x m] = —ian , [n x m] = —iom . (66) 



The set of vectors {n, m, rh} forms a basis in a 3- 
dimensional complex linear vector space. One can check 
that these vectors are the eigenvectors of L which cor- 
respond to the eigenvalues {1,0,2} respectively, and the 
operator L in this basis has the form 





1 








L = 



















2 



As expected, the operator L is degenerate and its matrix 
rank is 2. 

Using the expression it is easy to see that the 
solution of Eq. (|57|) is 

/ = am , (68) 

where a is an arbitrary complex quantity. Thus, one has 

(n,/ )=0. (69) 

The expressions (|39l) and imply that the vector n 
is parallel to the Poynting vector V in the geometrical 
optics approximation. 

It is easy to check that because of the degeneracy of the 
operator L, the following relation is valid for an arbitrary 
vector a: 

La = n(n, La) — ia[n x La] . (70) 

This relation shows that Eq. (|53[) has a solution only if 
its right hand side obeys the linear constraint 

curl/ = n(n, curl/) — ia[n x curl/] . (71) 

This constraint is exact. 



C. Ray trajectories 

The 'standard' eikonal Eq. J5SJ 

(VS - gf = 1 (72) 

does not depend on the state of polarization of the wave. 
In fact, one can see that it coincides with the eikonal 
equation for the massless scalar field <p, which obeys the 
equation Dip = 0. 

Equation (|72|) has the form of the Hamilton- Jacobi 
equation. To solve it one uses solutions of a dynamical 
system defined by the following Hamiltonian: 

U(x\ V i) = \{P - gf = \l lJ {p* - 9i)(Pj - 9j) , (73) 
and the canonical symplectic form 

3 

fi = H d P* A dx ' ■ ( 74 ) 

i=l 
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The corresponding Hamilton equations are 



dx 1 &H ... . 



9j) 



1 



~ ^liiPj ~ 9j)(pi ~ 9l) ■ 



(75) 



(76) 



The value of the Hamiltonian is constant on the ray tra- 
jectories. We choose 

1 

2 ' 



H 



(77) 



This condition determines the choice of the parameter 
along the trajectories. One has 



dx 

dC 



==r =1. 



(78) 



Hence i is the proper distance defined along an optical 
ray trajectory in the metric 7^. 

It is easy to show that the Hamilton equations (|75l) 
(1751 are equivalent to the following second order equa- 
tion: 

D 2 x 



dP 



dx 

- x curlg 



(79) 



where D/d£ is the covariant derivative defined with re- 
spect to the metric jij along the trajectory. 

Let = be an equation of a 2-dimensional sur- 

face r which defines the initial position of the wave front. 
Denote p^ = Vj$(a;). For each initial point x of V and 
the corresponding value p^ one can solve the Hamil- 
ton equations (f75 ]) ~ (j76)) . The solutions define ray tra- 
jectories which form a 3-dimensional subspace of the 6- 
dimensional phase space called a Lagrange manifold. Its 
characteristic property is that the symplectic form re- 
stricted to this subspace vanishes. The action S calcu- 
lated on these trajectories is a function of the final point 
x, 



S(x)= / (p,dx) 



Xq 
X 



Xq 



1 



dx 
9, -77 



dt. 



(80) 



In the 3-dimensional configuration space such 2- 
dimensional family of rays covers a 3-dimensional region. 
We assume that through each point x passes only one 
ray of this family, so that the expression (j5H)) defines the 
function S(x). This function is the required solution of 
the Hamilton- Jacobi equation (f?2"j) (see, e.g., [25U27I ] ) . 

Before we discuss properties of the amplitude / let us 
make two useful remarks. 

1. Equation ([?9"f is analogical to the equation of mo- 
tion of a nonrelativistic particle of the unit charge- 
to-mass ratio in a 3-dimensional space in the pres- 
ence of a magnetic field defined by the vector po- 
tential g. 



2. This equaton is nothing but a (3 + l)-dimensional 
form of the 4-dimcnsional geodesic equation for a 
null ray in a stationary spacetime 



d\ 2 



0, xf 



fox*), 



(81) 



where A is an affine parameter along the null ray, 
and D/dX is the covariant derivative defined along 
trajectory of the null ray. Equation (|79|) can also be 
derived from (3+l)-decomposition of the geodesic 
equation (f8Tj) for null rays (see, e.g., [57|). 



D. Transport equations 

To derive the transport equations for the complex am- 
plitude / we apply the operator curl to Eq. (|53l) . Using 
Eq. (IMl) we derive 



curl/ — jcrcurl[n x /] = ui crcurlcurl/ . 



(82) 



If one substitutes the expansion (1551) into this equation 
the leading order term gives the following relation: 



curl/ = itrcurl[rj. x / ] . 



(83) 



This transport equation defines the complex amplitude 
f . Similar equations can be written for the higher order 
coefficients /j>i's. 

Using the properties of the operator L we constructed a 
local basis {n, m, fh}. The vector m is defined up to the 
transformation m — ¥ exp(i^p)m, where tp is a function of 
x. Let us choose the basis on the surface of the initial 
position of the wave front T and define it along the ray 
trajectories by the condition that it is Fermi propagated 
along the tangent to the rays vector n. Denote w = 
V n n, then the Fermi derivative of a vector a along a 
unit vector field n is 



J„o, = V n a — (n, a)w + (w, a)n . 



(84) 



One has J-nn = 0, and the scalar product of any two 
Fermi propagated vectors is constant. Since (m, n) = 
at the initial surface, this orthogonality condition is valid 
along ray trajectories. Thus, the condition that m is 
Fermi propagated takes the form 



V n m = —(w, m)n . 



(85) 



This condition fixes gauge ambiguity in the choice of the 
complex basis {n,m,rh}. We call such a choice of the 
basic vectors fields canonical. 

Denote /o = (/oi/o) 1 ^ 2 ' then the complex polariza- 
tion vector in the canonical frame takes the form 



(86) 



To find the amplitude fo and the phase if we proceed as 
follows. Using Eq. (|25|) we present Eq. (|83|) in the form 



curl/ = ia [ndiv/ - / divn + (/ , V)n - V n / ] , 

(87) 
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where V n = (n, V). Equations (|27| and ((57)) give 

div/ = ia[(f , curl n) ~ (n, curl/ )] . (88) 
and Eqs. (gBJ) and ([Ml) give 
(/„, V)n = -[nx curl/ ] - [/„ x curln] - V n / . (89) 



Substituting the expressions (|88| and (|89|) into (f87[) and 
using the relation (|7ip we derive the derivative of / 
along n 

V n /o = g [iem(/ , curln) - [/ ,curln] - / divn] . 

Taking the scalar product of this expression and the com- 
plex conjugate f and using Eq. (1231) we derive 



V n / 2 + / 2 divn = div(/ 2 n) = . 



(91) 



This expression implies conservation of number of pho- 
tons in a stationary spacetime. 

Using Eqs. ([55]). (|9"0|) , and flU]) we obtain the following 
relation: 

imV n j3 = - (ian{m, curl n) — [m, curln]) — V n m. 

(92) 

Analogously, we can define derivation of the vector n 
along itself. Using Eqs. (|26]l and ([59I we have 

V n n = — [n x curln] = [n x curlg] . (93) 

Then, using Eqs. (J66J) , flgSJ), j92]), and ([93l) we derive 

(94) 

Solving this equation we derive 



VnV' = -(n,curlg) 



aip, 



1 

V = o 



(curlgr, da;) . 



(95) 



As above, the integral is evaluated along a ray trajectory 
connecting a point on T and a point x. As a result, 
solving the transport equation along the 2-dimensional 
family of rays determined by the Lagrange manifold one 
finds the phase ip(x). 

To summarize, the application of the 'standard' geo- 
metrical optics gives the following expression for the vec- 
tor field T a : 



dx 



9 = 9 + 7T curl 9- 



(96) 
(97) 

(98) 



where we restored the superscript a. 

To finish this section consider a superposition of the 
left and right circularly polarized light beams of the same 
amplitude 



This field describes linearly polarized beam of light. Us- 
ing the expression (|96l) we find 

T « V2fok e^ s , (100) 
fc = ei costp — e 2 sin 93 , (k 07 k ) = l. (101) 

Here <p is determined along each of the ray by the rela- 
tion (|95p . These relations show that the vector of the 
linear polarization fco rotates with respect to the Fermi 
propagated basis with the angular velocity 



d<P 1, , s 
w = -(curlg,n). 



(102) 



This is the well-known gravitational analogue of the Fara- 
day rotation (see, e.g., [H [H [H, H3|). 

Let us mention that applying similar steps as above 
to the higher-order terms in Eq. (|56[) one can calculate 
the post-geometrical optics corrections to the obtained 
solution. 



V. MODIFIED GEOMETRICAL OPTICS 
A. Modified eikonal equation and ray trajectories 

We already mentioned that in the 'standard' geomet- 
rical optics the eikonal equation does not depend on the 
helicity. As a result, in such a formalism photons of differ- 
ent helicity propagate identically. The spin of a photon 
experiences a Faraday rotation, but an expected back- 
reaction of the spin on the photon's trajectory is absent. 
We already discussed this general problem in Sec. II. 
The origin of this problem is that the 'standard' geo- 
metrical optics is valid locally. However, to desribe the 
spin effects for photon scattering one needs to modify 
this approach. For this purpose one needs to include the 
lowest order helicity corrections into the eikonal equa- 
tion. Such a modification does not change the local re- 
sults, and therefore is in agreement with the 'standard' 
geometrical optics locally. However, at large distances, 
due to the accumulation of these small corrections, the 
trajectories of photons of different helicity are modified 
considerably. 

Equation (|96[) gives us hints how to made the required 
modification. In this representation the equation for the 
change of the amplitude /o along the ray trajectory does 
not depend on helicity. Moreover, the equation for the 
Fermi transport does not depend on helicity cither. The 
only helicity-dependent term is the phase <p° ' . But it is 
naturally combined with the geometrical optical phase S 
[see the relations (|97l) - (|98)) ]. In a general case, there is an 
ambiguity in the interpretation of such a phase factor. It 
is evident that the form of Eq. (|9"6")l remains the same if 
one makes the following transformation: 



f^fe 



iip 



S - oj-^j , 



(103) 



(99) of the vector amplitude / and the phase S. 
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The idea of the modified geometrical optics is to fix 
this ambiguity by imposing the condition that the trans- 
port equations for the vector amplitude /, at least in 
the leading order, do not depend on helicity. This has 
a very non-trivial consequences. Namely, in order to re- 
produce the expression for S (x) in this gauge, one needs 
to modify the eikonal equation. As a result, the effec- 
tive Hamiltonian for light rays will be also modified by 
helicity-dependent corrections. The corresponding La- 
grange manifold for our problem will be also modified. At 
small distances this modification is small, but at large dis- 
tances the corresponding corrections become important. 
Since the effective Hamiltonian for light rays and the 
Largrange manifold depend on helicity, they contain in- 
formation about the helicity influence on photon trajecto- 
ries. As expected, the corresponding helicity-dependent 
terms are suppressed by the small factor a; -1 . 

After these general remarks let us describe the modi- 
fied geometrical optics equations in detail. Let us rewrite 
the main equation (|55|) by adding the term icj -1 [curlg X 
/] to its both sides. The form of this modification is sug- 
gested by the effective eikonal jSTJ), 



The modified Hamilton equations are equivalent to the 
following second order equation for photon trajectories 
[cf. Eq. Cm)]: 



D 2 x 
~dP~ 



dx 



x curlg 



2uu 



dx 



x curl curlg 



(112) 



This equation explicitly shows the effect of the helicity 
on the ray trajectories. Repeating the steps of Sec. IV 
for g replaced with g we derive the effective eikonal [cf. 
Eq. (HH)] 



S(x) = 



x P 



di + (g, dx) + — (curl g, dx) 



(113) 



The main difference with the previous result is that the 
integration in (|113[) is performed over the modified ray 
trajectories corresponding to the vector n. The last term 
in this equation depends on the helicity sign a. 



B. Modified transport equations 



Lf + ^~ [curlg x /] = -curl/ + J- [curlg x /] . (104) 
Zuj lu Zuj 

Accordingly, we define the following new linear operator 
L [cf. Eq. (HI]: 



Lf=f-ia[fixf}=0. 



(105) 



n = p -g, p = VS\ g = g+— curlg. (106) 
Then, Eq. (|104j) takes the form 



Let us now discuss modifications of the transport 
equations. To begin with, let us introduce the basis 
{n, rh, rh} such that the relations between the basis vec- 
tors are analogous to the relations ([65| and (l66|) . Analo- 
gously, we define [cf. Eq. (1551) ] 



f = f rhe lv , 



(114) 



so that 



Lf = -curl/ + — [curlg x /] 

lu Zuj 



(107) 



Substituting the expansion (|55l) into this equation and 
taking the geometrical optics limit uj^ 1 — > in accor- 
dance with the large distances asymptotic approxima- 
tion, i.e., we keep the term of the first order in uj^ 1 which 
enters the operator L, we derive [cf. Eq. §51 



£/o = o. 



One can see that the operator L has the same prop- 
erties as the operator L. In particular, we have [cf. Eq. 



(h, h) = 1 . 



(»,/„) =0. (115) 

From the relation ([70]) . which remains valid for the op- 
erator L and the vector n, and Eq. (|107[) we derive the 
following relation: 



curl/ = h(n, curl/) — ia[fi x curl/] — — [curlg x /] 

(108) ia 1 1 
-yn(curlg, [n x /]) + -(n, /)curlg - -(n, curlg)/. 

(116) 

Applying the operator curl to Eq. (|107[) and using Eq. 

(109) (I105P we derive 



where h is the eigenvector of L with the eigenvalue 1. 
According to the definition of the vector n, the effective 
eikonal equation reads [cf. Eq. (|72l) ] 



(V5-g) 2 = l. 



(110) 



This modified eikonal equation can be solved as above. 
It coincides with the Hamilton- Jacobi equation for the 
following modified Hamiltonian: 



1 



1 



H(x\ Pl ) = ^(p-g) 2 ^ ^{ Pi - ~ 9i ){ Pj - g 3 ) , (111) 



(7 2 

curl/ — zcrcurl[n x /] = — curl curl/ 4- — curl[curl<7 x /] . 

oj Zuj 

(117) 

Substituting the expansion (|55[) into this equation we de- 
rive the following equation corresponding to the geomet- 
rical optics approximation [cf. Eq. ([83j> ]: 



curl/ = icrcurl[n x / 



(118) 



Repeating the same steps as in Sec. IV and using Eqs. 
(|114p . (|115l) . and the relation (|116j) we derive [cf. Eq. 
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C. Example: Kerr spacetime 



Vfi/ = - io-h(f , curl n) - [/ ,curln] - / divn 



T 



«(/o> curl 5) ~ io "[/o) curl 0] + /o(™,curlg) 



(119) 



Using this expression, one can show that the following 
relation holds [cf. Eq. (j9l~1) ]: 



div(/ 2 n) = 0. 



(120) 



That is, conservation of number of photons in a station- 
ary spacetime holds for the modified ray trajectories. Us- 
ing the expressions (|114[) . (|119p . and (| 1201) we obtain the 
following relation [cf. Eq. (|92p]: 

irhV^if = — (ian(rh, curl n) — [rh,curln]) 



— [n(ra, curl g) — ia[rh, curlg] 

+ rh(n, curl g)] — \7 A rh . 



(121) 



Analogously, using Eqs. ((26]), (|106[) . and f|109|> we can 
define derivation of the vector n along itself [cf. Eq. 

H23M, 

V^n = — [n x curln] = [fi, x curlg] + ^—[™ x curlg] . 

(122) 

We use the basis {n,m,m}, which is Fermi propagated 
along the corresponding ray trajectories generated by h. 
Then, using Eqs. (|T2"T|) and (H"2"2"j) we derive [cf. Eq. (IM1) ] 



— (n,curl curlg) 
4w 



(123) 



Thus, in the geometrical optics limit uj^ 1 — > we have 
(p = const, which implies that polarization vector does 
not rotate with respect to the basis {h,m,fh}. 

To summarize, the expression for the field J- a in the 
modified geometrical optics approximation reads 



5(x) = 



rX 




/_ dx\ 


lx 


1 + 





9 = 9 + —curie/. 



(124) 
(125) 

(126) 



This result is similar to (|96p , but it contans an important 
difference. To define the Fermi propagated basis and the 
integration paths one must use the modified ray trajec- 
tories defined by the tangent vector h. 

Since a linearly polarized light is a superposition of 
the right and left circularly polarized beams, we conclude 
that the initially linearly polarized beam in a stationary 
gravitational field splits into two beams with different 
circular polarization. However, it is easy to check that 
before this effect becomes important the beam split is 
negligible, and the linear polarizion vector has the same 
Faraday rotation (|102l) . 



Let us emphasize that in the developed scheme of the 
modified geometrical optics we did not assume that the 
gravitational field is weak. One can expect that spin- 
optical effects might become significant for the light prop- 
agation close to rotating black holes, where the gravity 
is evidently very strong. We shall discuss this subject in 
detail in the future publication. Just to give an exam- 
ple, here we present the terms which control the helicity 
dependence of light rays in the Kerr geometry. 

The Kerr metric in the Boyer-Linquist coordinates 
x M = (t, r, 9, 4>) can be presented in the form ^ with 



h = (A - a sin 2 9) /£ , g l = - 



2aMr 



(127) 



di 2 = j^dx'dx 1 = -^-dr 2 + jd9 2 + AS ™ 9 d<t> 2 , (128) 

L J It I L It 



where 



2Mr 



(129) 



This metric represents a Kerr black hole of the mass M 
and the angular momentum J = aM, where < |a| < M. 
For this metric we derive 



(cmlgf = ^ — cos 95\ 

2aM(r 2 - a 2 cos 2 i 



X 3 

(curl curl gf = 



sin f 



(130) 
(131) 



Note that Eq. (|112[) for photon trajectories contains an 
extra term (|131[) which in the weak field approximation 
gives the angular separation of the right and left circu- 
larly polarized beams of the order of GJX/(c 3 L 3 ), where 
L ^> GM jc 2 is the distance from the Kerr black hole. 
This result is in agreement with the estimation obtained 
by Mashhoon [ill. 



VI. SUMMARY 

In this paper we studied the propagation of circularly 
polarized monochromatic electromagnetic waves in a sta- 
tionary spacetime. We presented the Maxwell equations 
in (3+l)-form corresponding to (3+l)-decomposition of 
a stationary spacetime. The Maxwell equations for circu- 
larly polarized monochromatic waves reduce to the set of 
equations for the complex amplitude of such waves. This 
set of equations can be solved if a solution to the cor- 
responding master equation for the complex amplitude 
is found. We first described the 'standard' procedure of 
finding solutions to the master equation in the geomet- 
rical optics approximation. We used the main result of 
this approximation (1961) in order to develop an improved 
modified scheme. This scheme locally gives the same 
results as the 'standard' geometrical optics, but for the 
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scattering problem it allows one to include the helicity 
corrections to the photon equation of motion. It is inter- 
esting that technically, in order to obtain these modified 
equations, it is sufficient to modify the gravitomagnetic 
vector potential g by the well-defined correction which 
has the form | era; curl g. Let us emphasize that the 
obtained results are valid beyond the weak field approx- 
imation for the gravitational field, and hence, they allow 
one to study the polarization effects in strong gravity. 
Using the proposed scheme one can calculate the post- 
geometrical optics corrections to the obtained solution. 



One can interprete the obtained results in a slightly dif- 
ferent manner. As a result of spin-orbit interaction, pho- 
tons with different circilar polarization move in a slightly 
modified metric. This modification depends on the fre- 
quency of a photon u>. Such 'renormalization' approach 
allows one to conclude that due to the scattering of cir- 
cularly polarized 'white' light, which contains different 
frequencies, one may observe a peculiar 'rainbow effect'. 
Another effect is that linearly polarized beam of light 
propagating in a stationary gravitational field splits into 
two components of the right and left circular polarization. 
These and other polarization effects might be important 
for propagation of polarized photons near the horizon of 
rapidly rotating black holes. 
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Appendix A: Calculation of determinant 

Suppose that a is some scalar quantity, J is a unit and 
A is an antisymmetric 3x3 matrices, respectively. Then 
one has 

det(al + A) = a 3 - |tr(A 2 ). (Al) 
The antisymmetric matrix A can be written in the form 

where e i2 3 = e 123 = 1 and a 1 = \e^ k A ]k . Using this 
representation one obtains 



Thus, 



tr(A 2 ) = -2a 2 = -2a ia l . 



det(a/ + A) = a 3 + aaia 1 . 
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